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Heating of carriers in an intrinsic graphene under dc electric field is considered taking into account 
the intraband energy relaxation due to acoustic phonon scattering and the interband generation- 
recombination transitions due to thermal radiation. The distribution of nonequilibrium carriers is 
obtained for the cases when the intercarrier scattering is unessential and when the carrier-carrier 
Coulomb scattering effectively establishes the quasiequilibrium distribution with the temperature 
and the density of carriers that are determined by the balance equations. Because of an interplay 
between weak energy relaxation and generation-recombination processes a very low threshold of 
nonlinear response takes place. The nonlinear current-voltage characteristics are calculated for 
the case of the momentum relaxation caused by the elastic scattering. Obtained current-voltage 
characteristics show low threshold of nonlinear behavior and appearance of the second ohmic region, 
for strong fields. 

PACS numbers: 73.50.Fq, 81.05.Uw 
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I. INTRODUCTION 



Active study of graphene, taking place in recent years, 
is stimulated both by unusual physical properties of this 
gapless and massless semiconductor (see discussion and 
references in the review [l[) and by the need to study 
the possibilities of applications in electronics, see discus- 
sion in [2]. For such applications study of a response 
of nonequilibrium carries is important and very much 
needed. Both experimental and theoretical studies of 
the nonequilibrium electron-hole pairs in graphene un- 
der interband photoexsitation have been already per- 
formed, see [|[ and [3,0] (and references therein), respec- 
tively. Experimental investigations of the carries heating 
in graphene by dc electric field were carried out in re- 
lation with demonstrations of graphene field-effect tran- 
sistor (see @ and references therein). Essential heating 
apparently had taken place in experiments on current- 
induced cleaning of graphene Q- In respect to theoret- 
ical treatments of the problem, it was investigated the 
rate of energy relaxation of nonequilibrium carries, both 
analytically [8| and numerically, [9| and also it was dis- 
cussed a derivation of kinetic equation for a strong field 
[Tfj| . However, to the best of our knowledge, the ba- 
sic problem of carries heating by dc electric field have 
not been considered so far and a systematic investigation 
of the non-linear current-voltage characteristics was not 
performed yet. 
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In this paper, we study a heating of electron-hole pairs 
in an intrinsic graphene under dc electric field. To de- 
scribe heating it is taken into account both increase of 
the energy of carries in the field, E, and the follow- 
ing relaxation processes: (1) damping of the momen- 
tum due to elastic scattering on structural disorder, that 
is the most fast process ensuring weak anisotropy of 
the distribution function [U, (2) intraband quasi-elastic 
energy relaxation on acoustic phonons, (3) generation- 
recombination processes due to interband transitions in- 
duced by thermal radiation, and (4) intercarrier scatter- 
ing due to Coulomb interaction. In present study we 
consider the heating within two limiting regimes: (I) if 
the Coulomb scattering is unessential, when it it is possi- 
ble to neglect by the intercarrier scattering, and (II) the 
Coulomb-controlled case. In the latter case the Coulomb 
scattering imposes the quasiequilibrium distribution on 
carries, with characteristics defined by balance equations 
for the density and the energy. In present study we find 
the distribution function of carries for these regimes and 
investigate it dependence on applied field and the tem- 
perature of thermostat formed by phonons and thermal 
radiation. In addition, we analyze current-voltage char- 
acteristics and the density of nonequilibrium carriers. 

The character of obtained nonlinear response is deter- 
mined by a few peculiarities of the considered model. 
First, the velocity of carries does not increase with the 
energy as it is equal vwp/p, where vw = 10 s cm/s is 
the characteristic velocity of charged neutrino-like par- 
ticles, p is the 2D momentum. Second, the interband 
transitions are effectively excited by thermal radiation 
(the matrix element of transition is oc Vyy) they not only 
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change a carrier concentration but also participate in the 
energy relaxation. As a result, the noncquilibrium dis- 
tribution is formed due to interplay between weak en- 
ergy relaxation and generation-recombination processes 
(i.e., between interactions with phonon and photon ther- 
mostats). Third, the rate of energy relaxation sharply in- 
creases as p grows, where weak electron-phonon interac- 
tion takes place for slow carriers. Therefore for scattering 
on short-range defects, when the relaxation of momen- 
tum also grows with the energy, it is realized a sub-linear 
current- voltage characteristic with essential nonlinearity 
for weak fields while as field grows nonlinearity becomes 
weak. If the momentum relaxation becomes ineffective 
for high energy carriers (because of a finite-range disor- 
der) a super-linear increase of current takes place. In the 
region of strong fields a linear dependence of current from 
the field is realized (the second ohmic region of current- 
voltage characteristic) where the effective conductivity is 
smaller than the conductivity of equilibrium carriers for 
low temperatures, T, and the former becomes larger than 
the latter as T grows. Because the intrinsic graphene has 
maximum resistance (with respect to gate voltage), here 
the Joule heating is least effective so maximal electric 
field can be applied. 

The paper is organized in the following way. The basic 
equations describing the heating of carriers in an intrinsic 
graphene under dc electric field are presented in Sec. II. 
In Sec. Ill we examine the symmetric distribution func- 
tion of nonequilibrium electron-hole pairs for the cases / 
and (II), see above. The current- voltage characteristics 
are analyzed in Sec. IV. The concluding remarks and 
discussion of the assumptions used are given in the last 
section. Appendix contains the evaluation of the quasi- 
classical kinetic equation for electron-hole pairs under a 
strong dc electric field. 



II. BASIC EQUATIONS 

Nonequilibrium electron-hole pairs in the intrinsic 
graphene are described by coinciding distribution func- 
tions fe.hp = f P as their energy spectra are symmetric 
and processes of scattering there are identical (see Ap- 
pendix) . Therefore instead of system of kinetic equations 
for fep and fhp it is possible to consider the single kinetic 
equation as 



eE 



dfp 
dp 



E-m/Ip}, 



(i) 



where E is a strong electric field and the classical form 
of the field term is substantiated in Appendix. Here 
p is the 2D momentum and Jj{f\p} is the collision 
integral for the jth scattering mechanism, with index 
j = D, LA, R, and C correspond to the static disorder 
(D), the acoustic phonon scattering (LA), the radiative- 
induced interband transitions (R), and the carrier-carrier 
scattering (C), respectively. Integrals of elastic scatter- 



ing there were considered earlier, see [TT], [L2J and refer- 
ences therein, Jla{J\p} and Jr{/|p} were evaluated in 
H, and Coulomb scattering was considered in [l3[. For 
the distribution function defined by Eq. (UJ), the concen- 
tration balance equation is given as 



72EM/|P} = 0, 



(2) 



because interband transitions are forbidden not only for 
elastic scattering but also for Coulomb scattering, due to 
symmetry of the energy spectrum (l4^ . and for phonon 
scattering, due to s <C vw', here s is the velocity of sound. 
In Eq. ([2]) the factor 4 takes into account spin and valley 
degeneracies, L 2 is the normalization area. 

Taking into account that electrons and holes equally 
contribute to the current density, I, one obtains 



p 



(3) 



where the asymmetric part of distribution function, A/ p , 
is separated from the symmetric one, f p , by using the 
relation / p = f p +Af p . Last contribution is small if usual 



condition |eE|Tp^ m ' 1 <C p is satisfied, [lfj Here Tp {m> is 
the momentum relaxation time and A/ p is given as 



.fi(r. 



A/ P = 



P 



dp 



(4) 



For the case of a random potential U x characterized 
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by the correlation function (U x U x i) = Ud exp{ — [(x — 
x')/Z c ] 2 } with the averaged energy Ud and the correla- 
tion length l c , one obtains 



1 



VdP 



B(m) 



-h (z 2 ) , (5) 



where the decreasing with pl c /h form-factor ^(z) is 
written through the modified Bessel function I\[x) and 
the characteristic velocity Vd is introduced as Vd — 

ttU dl 2 /^ 2 vw- The conductivity, a is determined ac- 
cording to the standard formula I = crE. Using the 
above-introduced relaxation time (5) and Eqs. (3, 4), 
one transforms the conductivity as follows 



2vw 



nfi Vd 



dp 



^{pic/h) 



dfp 
dp 



(6) 



where the averaging over angle is performed. 

Further, we turn to the averaging over angle (below 
we symbolize such averaging as a line over expression) of 
Eq. ([1]). Neglecting the weak contribution of A/ p in the 
right-hand side of (fT]), one obtains the kinetic equation 
for symmetric distribution f p in the following form 



eE 



dAfp 
dp 



(J) 
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Here Jj{f\p} = Jj{t\p} and summation is performed 
over j = LA, R, C because the elastic scattering does 
not affect the symmetric distribution due to the energy 
conservation. Performing the averaging over the angle, 
we transform the field contribution of Eq. ([7]) as follows 



m 



eE- 



dA/p = (eg) 2 d 
dp 2p dp 



df P 
dp 



ir ■ (8) 



As a result, using the Fokker-Planck form of Jla{/\p} 
and Jtin{f\p} obtained in 0], we arrive to kinetic equa- 
tion 



v p qe) d 
p 2 dp 



dfp 
dp 



Pt 



- fp) 



+^ {R) [N 2p / PST (1 - 2/p) - fp] + -hc{f\p} = 0, (9) 

where N x — (e x — is the Planck distribution and 
Pt = T/vw The rates of quasi-elastic energy relax- 

, can be 



ation, v>. 



fi(ge) 
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and of radiative transitions, 



MR) 



presented in the form [5j 



V 



(-)' 

\vw J 



VacP 



D 2 T 



AH p s V\yS 2 ' 



MR) 



VrP 

h ' 



he 



8vw 



(10) 



The characteristic velocities v ac and v r are introduced 
here in order to separate the linear momentum depen- 
dence and expressed through the deformation potential 
and the sheet density of graphene, D and p s , as well as 
the dielectric permittivity, e. Characteristic momentum, 
Pe, is defined in Eq. ^ by relation 



4 = fvw\ 2 jeEh) 2 
PE \ s ) v ac v d 



(11) 



so that pe oc \/E. The Coulomb scattering integral 
Jc{f\p} can be neglected in Eq. ([9]) for the case I 
therefore here nonequilibrium distribution f p is defined 
by nonlinear differential equation of the second order. 

For other approach, in the case II, the Coulomb scat- 
tering term in Eq. Q is dominant and it imposes 
quasiequilibrium distribution 



f p = {cxp[(v w p - fx)/T c ] + 1}' 



(12) 



with effective temperature of carriers T c and quasichem- 
ical potential /i. To determine T c and p we will use the 
concentration balance equation ([2|) and the energy bal- 
ance equation. The latter we obtain by summing of Eq. 
([7]) over p with the energy weight vwP as follows 



J2p[ J LAif\p) + J R (f\p)} = 0, (13) 



where the field term is expressed through the conductiv- 
ity © using the integration by parts (the factor 1/2 is 



quite understood as the total Joule heat aE 2 must to 
be divided equally between electrons and holes). There- 
fore, the two parameters of f p are determined by two 
transcendent equations ([2]) and (fTB"]) . Now, the nonlinear 
conductivity should be calculated under substitution f p 
into Eq. (6). 



III. NONEQUILIBRIUM DISTRIBUTION 

Below we consider the nonequilibrium distribution ob- 
tained from Eq. ((9]) with the Coulomb contribution omit- 
ted (case /) or from the balance equations © and (fT3")) 
(case 77) . Also we present the field and temperature de- 
pendencies of the nonequilibrium sheet concentration. 



A. Weak intercarrier scattering 

We start with consideration of the case I, when the 
carrier-carrier scattering is ineffective. Omitting Jc in 
Eq. ([9]) and introducing the dimensionless momentum 
x = pIpt and the parameter 77 = prlc/h, we obtain the 
nonlinear differential equation 



d 

d.r 



(pe/pt) 4 



2* (x T]) 
+Tx 2 



f +*/.<!-/.: 



1-2/, 

e 2x _ I 



-II =0. 



(14) 



Here the dimensionless parameter T = {vw/s) 2 v r /v ac 
determines the ratio of contributions from the thermal 
radiation and acoustic phonons, to the energy relaxation. 
As the boundary conditions for Eq. (fT4Tj we use that 
x (dfx/dx + f x )x^oc — 0, so that f x must decrease for 
x — > 00 sufficiently fast, and also the density balance 
equation ((2]). [15] If to integrate Eq. (fT4| over x from 
to 00, then for finite electric field besides Eq. ^ we 
obtain additional, proportional to (pe /PT) 4: (dfx/dx) x= o, 
term that, certainly, must be equal to zero. Thus, instead 
of Eq. @, the second boundary condition can be used 
in the following form {df x /dx) x ^o = 0. 

The numerical solution of Eq. (fT4|) is performed below 
with the use of the finite difference method and the iter- 
ations over non-linear contributions, see [l7[. Distribu- 
tion functions obtained are presented in Figs. 1 and 2 as 
functions of the energy for different temperatures, electric 
fields, and l c . For low temperature region, two regimes 
of the modification of a distribution emerge with the in- 
crease of field: weak electric fields, (pe /pt) 4 /(4r) <C 1, 
and strong electric fields, (p e I 'pt) a I '(4T) ^ 1- 

As it is seen from Figs, la and 2 in weak fields there is 
essential suppression of the distribution in the region of 
small p < pt (as slow carriers are heated very effectively) 
and for large p > pt distribution remains almost equilib- 
rium one. Further, as the field becomes strong, Figs. 1 
and 2 show that the carriers are spread over a wide region 
of energies such that, e.g., at T=4.2 K for E >5 V/cm the 
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FIG. 1: (Color online) Determined from Eq. (|14[) distribution 
functions, at T=4.2 K and l c = 10 nm: (a) for weak electric 
fields E =0.1 mV/cm (1), 1 mV/cm (2), 10 mV/cm (3), and 
30 mV/cm (4); (b) for strong electric fields E =0.1 V/cm (5), 
0.3 V/cm (6), 1 V/cm (7), and 5 V/cm (8). Dotted curves 
correspond to the equilibrium distribution. 



tail of the distribution will reach the energies for which 
the spontaneous emission of the optical phonons begins 
(^90 meV). It is seen, as well in agreement with Figs. 1 
and 2, that with the increase of T the transition between 
the weak field and the strong field regimes occurs at a 
higher E oc T 2 . In Figs. 2a - 2c we plot the distribution 
functions for l c = lOnm and 20nm for different E and 
T. It is seen that with the increase of E the influence 
of a finite l c grows, by making the effect of electric field 
on the distribution function more pronounced for higher 
l c . However, for weak E the difference between the solid 
curve and it dashed counterpart is practically invisible. 
Notice, from the solid and, especially, the dashed curves 
4 in Fig. 2c it follows that here interaction with optical 
phonons can be important. 



B. Coulomb controlled distribution 

Further, we examine the case //, when quasiequilib- 
rium distribution Eq. (|12| is determined from the bal- 
ance equations, ((2|) and (fl3l) . Introducing dimensionless 
momentum, y — v\yp/T c , such that f y — [exp(y— /x/T c ) + 
1] j we rewrite the concentration balance equation ^ 
in the form 



dyy 



1-2/, 

e 2yTJT _ l 



f 
■I v 



0. 



(15) 



0.3 
0.2 
0.1 
0.0 

0.4 

0.2 

0.0 
0.4 

0.2 

0.0 










(b) 






^3 


4 









10 20 30 40 50 60 









(c) 








< i < i < i < i < i < 



20 40 60 80 100 120 

u w p (meV) 



FIG. 2: (Color online) The same as in Fig. 1 calculated: (a) 
for T=77 K and E =0.14 V/cm (1), 1.4 V/cm (2), 4.4 V/cm 
(3), and 14 V/cm (4); (b) for T=150 K and E =74 mV/cm 
(1), 0.74 V/cm (2), 2.35 V/cm (3), and 7.4 V/cm (4);(c) for 
T=300 K and E =1 V/cm (1), 3 V/cm (2), 10 V/cm (3), and 
30 V/cm (4). The solid and the dashed curves correspond to 
l c = 10 nm and 20 nm, respectively. Dotted curves correspond 
to equilibrium distributions. 



which gives the relation between /x/T and dimensionless 
temperature T c /T (this relation is not dependent explic- 
itly on the field). Now the energy balance equation (Qji 
is given as 



T 



T 







1 - Vv 



0. 



(16) 



where field contribution Qe is transformed using Eq. ([6]) 

4 



Qe = 



PE 



T c /v w 



fy=0 + 



dyfy^ivcy) 



(17) 



Here r\ c = T c l c /hvw and it is introduced the func- 
tion $(z) = — ty'(z)/^(z) 2 and taken into account that 
*(0) = 1/2. 

In Figs. 3a and 3b calculated from Eqs. (fl~5]) . fTB]) 
dimensionless effective temperature T c /T and maximal 
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value of the distribution function f p —o = [exp(— /i/T c ) + 
1] j that determines quasichemical potential, are shown 
as function of E for different T. From Fig. 3a it is seen 
that T c /T grows faster with the increase of E for larger l c ; 
in addition, the grows of dimcnsionlcss effective temper- 
ature becomes faster for smaller T. Figure 3b shows that 
the characteristic value of distribution, f p= o, decreases 
faster with E at smaller T, 
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FIG. 3: (Color online) Determined from the balance Eqs. 
(|15|) . (|16|) dimensionless effective temperature T c /T (a) and 
maximal distribution f p =o (b) versus electric field for T =20 
K (1), 77 K (2), 150 K (3), and 300 K (4). Solid and dashed 
curves are correspondent to l c = 10 nm and 20 nm, respec- 
tively. 



C. Hot carrier concentration 

Using the solutions of Eq. I|14p or of the balance 
Eqs. (fl"5|) - (fl"7|) . we calculate below the sheet concentra- 
tion which is determined as follows 
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FIG. 4: (Color online) Normalized carrier concentration, 
n/riT, versus electric field: (a) at T =4.2 K (1), 20 K (2), 
77 K (3) and (b) at T =150 K (4), 300 K (5). The solid and 
the dashed curves are calculated from Eq. (|14p for l c =10 
nm and 20 nm, respectively. The dotted and the dot-dashed 
curves are calculated from Eqs. (|15p -(|17 p for l c =10 nm and 
20 nm, respectively. 



K, 150 K, and 300 K one obtains nr = 1.6 x 10 7 cm 2 , 
3.6 x 10 8 cm" 2 , 5.4 x 10 9 cto" 2 , 2.0 x 10 10 cto- 2 , and 
8.1 x 10 10 cm~ 2 , respectively]. From Fig. 4 it is seen 
that: i) at low temperatures (and not too large E) the 
relative increase of the density (n — nrj/iij- is approxi- 
mately oc E 1 / 2 as due to small characteristic momentum, 
p, effect of a finite l c is negligible; ii) at higher temper- 
atures (or very large E) (n — n T )/n T w A + BE, i.e., 
it is linear as due to large p effect of a finite l c becomes 
essential. Notice, that the density grows becomes faster 
for larger l c . 



IV. CURRENT- VOLTAGE CHARACTERISTICS 



n = 



dppfp 



(18) 



Here in the right hand side f p is replaced by f p for the 
case 77. As a result, the field dependence of n is deter- 
mined by a competition of the effective temperature T c 
grows and the maximal distribution / p =o decrease, see 
Figs. 3a and 3b: so that the concentration will grow 
slowly with the field. 

In Fig. 4 we plot the dimensionless carrier con- 
centrations n/riT as function of E. Here the equilib- 
rium concentration readily follows from Eq. (fT5)) as 
n T = n(T/hv w ) 2 /6 [notice, at T = 4.2 K, 20 K, 77 



Using obtained in Sees. IIIA and IIIB nonequilibrium 
distribution functions we calculate here nonlinear con- 
ductivity introduced by Eq. ^ and analyze modifica- 
tions of the current voltage characteristics due to tem- 
perature and l c variations. Performing the integration in 
Eq. ([6|) by parts one arrives to 



a = a 



lr 



dpf p <&(pl c /h) 



(19) 



where cto = (2vw/vd)e 2 /irh is the characteristic conduc- 
tivity. For the case of short-range scattering, pl c /K <C 1 
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FIG. 5: Dimensionless conductivity, cr/oo, versus field E for 
lc =10 nm at different temperatures: T — 4.2 K (a), 20 K 
(b), 77 K (c), 150 K (d), and 300 K (e). The solid and the 
dashed curves are calculated from Eqs. (|14f) . (I19p and Eqs. 
(|15l) -(|17 [ ). (|19|l . respectively. Inset in panel (a) shows low- field 
dependencies. 



(p is the characteristic momentum of hot carriers) con- 
ductivity can be expressed through the distribution func- 
tion of low energy carriers as: a ~ 2ero/ P =o- As f p =o = 
1/2 for E — ► 0, it follows that <to is the linear (in the ab- 
sence of heating) conductivity in short-range scattering 
limit. In Fig. 5 we plot ct/cto versus field E for l c =\Q 
nm. Notice, in Figs. 5a-5c, except the inset of Fig. 5a, 
the same field region is used; it is different from the one 
used in Figs. 5d-5e. It is seen that both approaches give 
similar dependences even for very strong field, for given 
T and l c . The dependencies ct/cto vs E for l c ~2Q nm 
are plotted in Fig. 6, where conductivity grows faster for 
higher temperatures. 

Next, we turn to consideration of the current-voltage 
characteristics (i.e., the current density / = aE versus 
field E). In Fig. 7 we plot the current density, / = aE, 
versus field E for Z c =10 nm at different temperatures. 
Similar dependencies for l c —2Q nm are shown in Fig. 8. 
It is seen that current-voltage characteristics can be es- 
sentially nonlinear already for relatively small E and / 
starting from 1CP 5 A/cm in the low-temperature region. 
Point out that for larger E and I (if E > 1 V/cm at 
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FIG. 6: The same as in Fig. 5 for l c =20 nm. 



low temperatures and E > 10 V/cm at room temper- 
ature) Figs. 7 and 8 will present dependences that are 
close to linear ones, A + BE. In addition, Figs. 7 and 8 
show that from both approaches (cases / and II) similar 
current-voltage characteristics are obtained, for given T 
and l c . It is also seen that current- voltage characteristics 
can manifest as sublinear so superlinear form in wide E 
regions depending on value of T, Z c ; these modifications 
of the form are essentially related with the effects due to 
finite L. 



V. CONCLUSIONS 

To summarize, we developed the theory of the carries 
heating in intrinsic graphene under strong dc electric field 
for the cases when intercarrier scattering is negligible or 
dominant. It is found that the deviation from the equilib- 
rium distribution starts from very low electric fields (from 
0.3 V/cm for the room temperature and from 10 fiV/cm 
for the liquid helium temperature; in particular, for a fi- 
nite fields dfp/dp tends to zero for p <C pe ^ Pt, (while 
for equilibrium distribution dj^ q jdp ~ —1/Apx has the 
finite constant value) due to ineffective energy relaxation 
of slow carriers. Since with the increase of heating this 
relaxation sharply grows and the recombination of car- 
riers is dominant over the thermogeneration within the 




E (V/cm) 

FIG. 7: Current-voltage characteristics for l c =10 nm at dif- 
ferent temperatures: T = 4.2 K (a), 20 K (b), 77 K (c), 150 
K (d), and 300 K (e). The solid and the dashed curves are 
calculated from Eqs. (fll)l. flj} and Eqs. p5 ]l -l[T7 ) . ([19]|. 
respectively. 

high energy region, it follows that the current-voltage 
characteristics should be sublinear for the short-range 
scattering case. The superlinear dependency that is re- 
alized in strong fields for scattering on finite-range disor- 
der. Therefore, for intrinsic graphene we obtain unusual 
combination of the low threshold of nonlinearity and ap- 
pearance of the second ohmic region for strong fields. 

Let us discuss the main assumptions made. Here it is 
studied intrinsic graphene which has at most resistivity 
so the effect of heating manifests itself the most. Doped 
materials require of separate investigation, moreover non- 
linearity in these materials must be weaker. We have re- 
stricted ourselves by study only of limiting cases of the 
absence of intercarrier scattering or dominating Coulomb 
scattering that imposes quasiequilibrium distribution. 
Point out, that the field dependences of concentration 
and the current-voltage characteristics there are in good 
agreement (the strict solution of the problem should be 
somewhere between the solutions obtained), i.e., they 
have weak sensitivity to the details of distribution func- 
tion. As the main generation-recombination mechanism 
it is assumed radiative-induced interband transitions be- 
cause Auger-processes are forbidden due to the symme- 
try of electron- hole states flij . Possible contribution of 




FIG. 8: The same as in Fig. 7 for l c =20 nm. 



other generation-recombination mechanisms (e.g., disor- 
der induced transitions, caused by acoustic phonons or 
carrier-carrier scattering) demands additional investiga- 
tion. At last, we consider interaction of carriers with 
equilibrium thermostat formed by acoustic phonons and 
radiation (compare with Ref. 18). The problem of heat 
removal is beyond of the scope of this paper. We specu- 
late now that heat removal is sufficiently effective [7| • In 
addition, it is possible to use short pulses (all times of 
relaxation there are < 1 /its), when the thermostat will 
not be overheated. These limitations are related with the 
lack of data on graphene, however, the qualitative picture 
will not essentially modify after taking into account with 
better precision of the parameters and the mechanisms 
of relaxation. 

The rest of assumptions there are standard. To de- 
scribe the momentum relaxation it is taken into ac- 
count only the statical disorder scattering using the phe- 
nomenological model of Ref. 11 (despite that, it corre- 
sponds well the experimental data, the microscopic mech- 
anisms of scattering are still unclear [19( ) and small con- 
tribution of acoustic phonons is discarded. Using of the 
quasielastic approximation for describing of the energy 
relaxation is justified by condition v\y ^ s, that ensures 
the small energy transference for the scattering process. 
As the energy of optical phonons is large, their contri- 
bution can be neglected even in considered strong fields. 



8 



Also it is unimportant in considered range of the ener- 
gies an anisotropy of both electron and phonon spectra. 
Interaction with thermic radiation can be described by 
taking into account only of direct interband transitions 
(the Drude absorption is small). Therefore, applied ap- 
proximations accurately describe the heating mechanism 
in graphene and lead to correct quantitative estimation 
of current-voltage characteristics. 

In conclusion, prospects of many device applications, 
in particular, the characteristics of field-effect transis- 
tors or effectiveness of graphene interconnections, are 
substantially dependent on the carries heating. In ad- 
dition, study of hot electrons (holes) gives information 
on electron-phonon coupling and about the mechanisms 
of recombination, i.e., makes possible verification of the 
relaxation mechanisms. Therefore, we believe that the 
obtained results will stimulate further experimental and 
theoretical study (as well as numerical modeling) of hot 
carriers in graphene. 



APPENDIX: KINETIC EQUATION 

Below we evaluate the system of quasiclassical kinetic 
equations governing the nonequilibrium carrier distribu- 
tions in graphene placed in a strong electric field E. Such 
a consideration is analogous to the approach used for the 
bulk narrow-gap semiconductors, see [2(|) and the recent 
papers [ici| which develop a similar approach for the case 
of graphene. But the case of the bipolar electron-hole 
plasma is not analyzed in these papers. 

We start from the single-particle density matrix which 
is governed by the quantum kinetic equation [la ]: 



iw 



eE • x, p 



J, 



coll ■ 



(A.I) 



The collision integral J co ii represents studied above scat- 
tering mechanisms and h\y is the 2x2 Weyl-Wallash 
Hamiltonian, [2l| describes states nearby the band cross- 
point. Here we use the linear dispersion laws ei p = IvwP, 
with I corresponding to the conduction (l=+l) or the va- 
lence (i=-l) bands, the eigenvectors \lp) are defined from 
the eigenstate problem: 

h w \lp) = £i P \lp), hw =vw(cr -p), (A. 2) 



lp> = 7i 



o»</> 



- lp> = 75 



where <fi is the p-plane polar angle. Distribution function 
over Ip states Fi p — (/p|p|ip) and the non-diagonal part 

of density matrix F p = (lp\p\ — lp) = (— lp|/3|lp)* are 
obtained from the system of kinetic equations: 



eE- 



OF, 



dp 



{f p x_« iI (p)-x,,_,( p )f;} 

= (lp\J c (p)\lp) (A.3) 
0F o 



-2v w pF p + eE- 
n op 



E-Xi,_i(p) 



x (F lp - F_ lp ) + -E • {X_ 1 ,_ 1 (p) - X M (p)} F p 

= (+lp|J c (/3)|-lp) (A.4) 

Here the interband matrix element of coordinate, 
Xj j'(p) = (Zp|x|Zp), is calculated on wave functions of 
the momentum representation (A. 2) and has the value of 
the order of h/p, where p is the characteristic momentum 
of the hot carries. 

In order to show the smallncss of nondiagonal com- 
ponents of the distribution function, the estimation 
(lp | J c {p)\ lp) ~ Fi p /r m is also used, where r m is the 
characteristic time of the momentum relaxation (the 
smallest characteristic time of the problem). For typical 
conditions of the applicability of quasiclassical descrip- 
tion [H 



2vwP > — , eEr c <C p 



(A.5) 



from Eq. (A.4) we have estimation F p /Fi p ~ 
h/(2v w pT c )- 

Further, neglecting by small nondiagonal contributions 
we arrive to the usual system of kinetic equations: 



eE- 



dF, 



ip 



dp 



3 



(A.6) 



which describes the distribution of hot carriers over the 
bands I = ±1 and the 2D momentum p. The scattering 
integrals, Jj{F\lp}, are obtained in Ref. [EG3- Here 
I corresponds to conduction (1—+1) or valence (1—-1) 
band, p is the 2D momentum, Jj{F\lp} is the collision in- 
tegral for the jth scattering mechanism, [j = D, LA, R, 
and C, see discussion after Eq.(l)]. It is convenient to 
make transition to the electron-hole representation in- 
troducing the electron (e) and hole (h) distribution func- 
tions, f ep and / ep , according to the replacements [22j : 



F+i,p ft 



ep j 



1 - F_i 



p' 



fi 



hp 



(A.7) 



and to rewrite the collision integrals in (A.6) through 

fe,hp' 

Finally, using the substitution (A.7) and the velocity 
operator v = i[hyy,x]/h = vy/cr, we obtain the velocity 
of lp state: (Zp|v|Zp) = lvwp/p = v; p . Taking into 
account four-fold degeneracy of states in graphene due 
to spin and valley degrees of freedom, one obtains the 
current density 



4e x - 
i P 



Aevw \ ^ P 



p 



(fep + hp) ■ (A.8) 



In addition, the intrinsic material the electron density is 
equal to the hole one (the neutrality condition): 



T5 ^2 (f £ P _ ^p) - 0. 
p 



(A.9) 



Thus, the symmetric scattering for the c- and the v- 
bands Eqs. (A6), (A8) and (A9) preserve their form 
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when / ep is replaced by fh P - Thus, the electron and / ep = /; lp = / p , so that the kinetic equation and the 
hole distributions in the intrinsic material are identical, current density take forms (1) and (3), respectively. 



[1] A.H. Castro Neto, F. Guinea, N.M.R. Peres, K.S. 
Novoselov, and A.K. Geim, Rev. Mod. Phys. 81, 109 
(2009). 

[2] P. Avouris, Z. Chen, and V. Perebeinos, Nature Nan- 

otechnology, 2, 605 (2007); P. Avouris, J. Chen, Materi- 

alstoday, 9 46 (2006). 
[3] D. Sun, Z.-K. Wu, C. Divin, X. Li, C. Berger, W. A. de 

Heer, P. N. First, and T. B. Norris, Phys. Rev. Lett. 101, 

157402 (2008). 

[4] S. Butscher, F. Milde, M. Hirtschulz, E. Malic, and A. 

Knorr, Appl. Phys. Lett. 91, 203103 (2007). 
[5] F.T. Vasko and V. Ryzhii, Phys. Rev. B 77, 195433 

(2008). 

[6] M.C. Lemme, T.J. Echtermeyer, M. Baus, and H. Kurz, 
IEEE Electron Device Lett. 28 282 (2007); Y. Q. Wu, 
P.D. Ye, M.A. Capano, Y. Xuan, Y. Sui, M. Qi, J.A. 
Cooper, T. Shen, D. Pandey, G. Prakash, and R. Reifen- 
berger, Appl. Phys. Lett. 92, 092192 (2008); X. Wang, 
Y. Ouyang, X. Li, H. Wang, J. Guo, and H. Dai, 
larXiv:0803. 34641 

[7] J. Moser, A. Barreiro, and A. Bachtold, Appl. Phys. Lett. 
91, 163513 (2007). 

[8] E. H. Hwang, B.Y.-K. Hu, and S. Das Sarma, Phys. 
Rev. B 76, 115 434 (2007); W.-K. Tse, S. Das Sarma 
larXiv:0812.1008l , 

[9] A. Akturkaand N. Goldsman, J. Appl. Phys. 103, 053702 
(2008). 

[10] M. Auslender and M. I. Katsnelson, Phys. Rev. B 76, 
235425 (2007); S.V. Syzranov, M.V. Fistul, and K.B. Efe- 
tov, Phys. Rev. B 78, 045407 (2008). 

[11] F.T. Vasko and V. Ryzhii, Phys. Rev. B 76, 233404 
(2007). 



[12] T. Stauber, N. M. R. Peres, and F. Guinea, Phys. Rev. 

B 76, 205423 (2007); N. M. R. Peres, J. M. B. Lopes dos 

Santos, and T. Stauber, Phys. Rev. B 76, 073412 (2007); 

F. Guinea. Journ. Low Temp. Phys. 153, 359 (2008). 
[13] L. Fritz, J. Schmalian, M. Mueller, and S. Sachdev, Phys. 

Rev. B 78, 085416 (2008); M. Mueller, L. Fritz, and S. 

Sachdev, ibid. 78, 115406 (2008). 

[14] M.S. Foster and I.L. Aleiner, larXiv:0810.4342l 

[15] F. T. Vasko and O. E. Raichev, Quantum Kinetic Theory 

and Applications (Springer, N.Y. 2005). 
[16] The relation A p +p(dA p /dp)/2 = [d(p 2 A p )/dp]/2p is used 

in Eq. (8). 

[17] D. Potter, Computational Physics (J. Wiley, London, 
1973). 

[18] P.N. Romanets, F.T. Vasko, and M.V. Strikha, Phys. 
Rev. B 79, 033406 (2009). 

[19] T.M. Mohiuddin, L.A. Ponomarenko, R. Yang, S.M. 
Morozov, AA. Zhukov, F. Schedin, E.W. Hill, 
K.S. Novoselov, M.I. Katsnelson, and A. K. Geim, 
larXiv:0809.TT62] S. Ada m, P.W. Brouwer, and S. Das 
Sarma, larXiv:0811.0609l 

[20] I.M. Dykman, V.M. Rosenbaum, and F.T. Vasko, Phys. 
Stat. Sol. (b) 88, 385 (1978) 

[21] E.M. Lifshitz, L.P. Pitaevskii, and V.B. Berestetskii, 
Quantum Electrodynamics, (Butterworth-Heinemann, 
Oxford 1982); P.R. Wallace, Phys. Rev. 71, 622 (1947). 

[22] A.I. Anselm, Introduction to Semiconductor Theory 
(Prentice-Hall, Englewood Cliffs, NJ, 1981); since there is 
no magnetic field under consideration, we do not reverse 
a momentum direction and the sign of charge (e < 0). 



